Measurement of the nonlinear Meissner effect in superconducting Nb films using a 
resonant microwave cavity: A probe of unconventional pairing symmetries 
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We report observation of the nonlinear Meissner effect (NLME) in Nb films by measuring the 
resonance frequency of a planar superconducting cavity as a function of the magnitude and the 
orientation of a parallel magnetic field. Use of low power rf probing in films thinner than the London 
penetration depth, significantly increases the field for the vortex penetration onset and enables 
NLME detection under true equilibrium conditions. The data agree very well with calculations based 
on the Usadel equations. We propose to use NLME angular spectroscopy to probe unconventional 
pairing symmetries in superconductors. 

PACS numbers: 74.25.-q, 74.25. Ha, 74.25.Op, 74.78.Na 



Meissner effect is one of the fundamental manifesta- 
tions of the macroscopic phase coherence of a super- 
conducting state. Meissner screening current density 
J = -en s v s induced by a weak magnetic field is pro- 
portional to the velocity v s of the condensate. At higher 
fields, the superfiuid density n s becomes dependent on 
v s due to pairbreaking effects, resulting in the nonlin- 
ear Meissner effect (NLME) [1-5] . NLME has attracted 
much attention since it probes unconventional pairing 
symmetries of moving condensates, for example, the d- 
wave pairing in cuprates or multiband superconductivity 
in pnictides. For a single band isotropic type-II super- 
conductor, the NLME at v s <C h/m^o is described by, 

J = -7^%(l-<o 2 Q 2 ), (1) 

where Ao is the London penetration depth, £0 is the co- 
herence length, Q = v s /mh = V9 + 27rA/0o, m is the 
quasiparticle mass, 9 is the phase of the order parameter, 
A is the vector potential, 0o is the flux quantum. 

NLME can manifest itself in a variety of different be- 
haviors, which reveal the underlying pairing symmetry. 
For instance, Eq. (1) describes a clean d-wave supercon- 
ductor at high temperatures fcsT > pfv s or a d-wave 
superconductor with impurities, where pf is the Fermi 
momentum [2-4], while in the clean limit at fc^T < Pfv s , 
the nonlinear term in Eq. (1) takes the singular form 
— [2], In the s-wave clean limit, the NLME is 

absent at T <C T c where a cx exp(— A/T) [6], but in the 
dirty limit a ~ 1 even for T — > 0. In multiband supercon- 
ductors NLME can probe the onset of nonlinearity due 
to the appearance of interband phase textures suggested 
for MgB2[7] or the line nodes and interband sign change 
in the order parameter or mixed s-d pairing symmetries, 
which have been discussed recently for iron pnictides [8] . 

To date, experiments aiming to observe the NLME 
in high-T c cuprates and other extreme type-II supercon- 
ductors have been inconclusive [9-11] mostly because of 



a very small field region of the Meissner state. Since 
NLME becomes essential in fields H of the order of 
the thermodynamic critical field H c = 0o/2v / 2VA t o^o£o 
[1] , penetration of vortices above the lower critical field 
HoHd = (</>o/47rAg)(ln k + 0.5) <C ^qH c imposes the 
strong restriction H < H c \, which reduces the nonlin- 
ear correction in Eq. (1) to ~ (£T cl /if c ) 2 ~ k -2 -c 1, 
where k = Ao/£o is the Ginzburg-Landau parameter. 
Yet even small NLME terms in Eq. (1) result in inter- 
modulation effects [3] under strong ac fields, as observed 
in YBa2Cu307_2: [12]. However, it remains unclear to 
what extent the intermodulation probes the true equilib- 
rium NLME, not masked by the nonequilibrium kinet- 
ics of gapless nodal quasiparticles excited by sufficiently 
strong ac field [13]. These features may contribute to 
some insufficiencies in experimental attempts to observe 
NLME both in magnetization [10, 11] and intermodula- 
tion [12] experiments. 

In this Rapid Communications we report the obser- 
vation of NLME, using a method which resolves the 
problems of vortex penetration and nonequilibrium ef- 
fects. This opens up the opportunity of probing NLME in 
any type-II superconductor under true equilibrium con- 
ditions. In the method illustrated by Fig. 1 the reso- 
nance frequency / = 1/2tt \[LC of a thin film planar 
cavity is measured as a function of a parallel dc mag- 
netic field Hp. Here C is the strip-to-ground capacitance, 
L = L g + Lfc is the total inductance, L g is the geomet- 
rical inductance, and Lk(H p ) is the field-dependent ki- 
netic inductance of the superconducting condensate mea- 
sured by this technique. For a film of width w, length s 
and thickness d -C Ao, we have Lk = ^asGX^/wd (Ref. 
[14]) where G ~ 1 is a geometrical factor [15]. The 
field- induced increase of X(H) yields the frequency shift 
Sf(H p ) = -SL k f/2L or 

51 If = -fi sG[X 2 (H p ) - \(0) 2 }/2Ldw. (2) 

The variation of is detected here using a very low 
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rf level inside the cavity (estimated to be ~ 10 pW 
at resonance), which makes it possible to measure the 
field-induced shift with high precision and sensitivity 
~ 1(T 3 - 1(T 4 in the linear response mode. Another 
key feature of this approach is the use of thin films of 
thickness d < Ao for which p,oH c i = (2(j) /wd 2 ) ln(d/£o) 
can be much higher than the bulk [J>oH c i [16]. As 
a result, NLME becomes much more pronounced as 
stronger fields can be applied without masking NLME 
by penetration of vortices. The NLME enhancement 
factor in a thin film can be evaluated from Q(z) — 
2-7r/ioAo J(z)/(f>Q induced by the Meissner current density, 
J(z) — H sinh(z/Ao)/Ao cosh(d/2Ao), which gives the 
maximum Q at the surface Q(^ lm = ird^H^ 1 " 1 /(f) Q for a 
thin film with d < A , and Q% lk = 2nX p, H h c ^ lk /<j> for a 
thick film d ^> Ao at their respective lower critical fields. 
Hence, the enhancement factor r = (Qm lm /Qm lk ) 2 = 
(dH f c l lrn /2\ H% lk ) 2 = [4A ln(d/£ )Mln K + 0.5)] 2 in- 
creases significantly as the film thickness decreases (for 
a YBCO film with d = 50 nm, A = 200 nm and 
£ = 2 nm, we obtain r ~ 10 2 ). The NLME con- 
tribution in Eq. 1 can be regarded as an effective in- 
crease of A 2 (H ) from the zero-field A 2 (0) to the maximum 
X 2 {H cl ) = A 2 (0)[l + a(4£ln(d/£)/d) 2 ]. Thus, NLME in a 
thin film is no longer cut off by low H c \, which can make 
the singular Yip-Sauls contribution [2] in d-wave super- 
conductors at ksT -C v s pf much more pronounced. 

Rotating the field in the plane of the film gives rise 
to an orientational dependence of X(H). This effect is 
illustrated by Fig. 1(b): the field H p inclined by the 
angle cp relative to the strip axis induces Meissner cur- 
rents J x and J y of which only J x along the strip cou- 
ples linearly with the weak rf current SJ U <C J x - For 
a wide film w 3> Ag/d, we have Q x — q(y,t) + pzsmtp, 
Q y = pz cos cp, p = 27T^o-ff p /0o and q(y,t) is due to the 
rf current along the film. Linearizing Eq. (1) in q and 
averaging over the film thickness, we calculate the weak 
rf current 6I U = —<pod J q{y)dy/2i:p J oX 2 {H) ) which de- 
fines the field-dependent X(H), 
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The constant a will be calculated below for the s-wave 
dirty limit. Here we emphasize two points: (1) rf currents 
excited in the stripline only probe the NLME caused by 
the strong dc Meissner screening currents uniform in the 
plane of the film; (2) the NLME correction is quadratic 
in d and in H p , and exhibits a cos2ip dependence, so 
that the NLME is minimum for the in-plane field applied 
along the strip axis. 

We demonstrate the angular dependence of NLME in 
a Nb film with k ~ 25. The film has a length s — 3 mm, 
width w = 100 /Km, thickness d = 65 nm, T c — 7.25 K, 
and the resistivity p n (T c ) — 23.2 /if2cm. The reduced T c 
is characteristic of Nb dirty thin films [17]. The length of 
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FIG. 1: (a) Cavity sketch as defined by the two RF ports 
(gaps) separated by half- wavelength \ mw /2 = 3 mm in an on- 
chip 50 SI coplanar waveguide. The unit vectors ftj.,2,3 show 
the coil axes, (b) Geometry of a thin film strip in a parallel 
field Hp, which produces the Meissner screening current loops 
depicted by dashed lines. 



the cavity corresponds to the half-wavelength resonant 
mode at ~20 GHz. The experimental setup is based 
on a heterodyne detector sensing microwaves transmit- 
ted through the two ports of the on-chip superconduct- 
ing cavity, as illustrated by Fig. 1(a) (similar microwave 
techniques have been implemented to detect X-rays in 
astronomy [18] or quantum states of superconducting 
qubits [19]). Measurements were performed at 80 mK 
in a Leiden Cryogenics dilution refrigerator. 

The challenge with the observation of NLME in a film 
is to eliminate the significant contribution from perpen- 
dicular vortices caused by the field misalignment with 
the film surface. This was achieved by aligning the film 
with the magnetic fields of three coils which produce a 
superposition of the orthogonal fields Hi oriented along 
the unit vectors ni.2,3 depicted in Fig. 1(a). To coun- 
teract the film misalignment by a proper field rotation, 
we performed measurements of the frequency shift Sf as 
a function of the nearly perpendicular field H 3 for dif- 
ferent in-plane field components (H\,H-i). The observed 
linear dependence of 5f{H 3 ) shown in Fig. 2 is consis- 
tent with the conventional behavior of A 2 (H ) determined 
by rf dynamics of pinned perpendicular vortices at low 
fields for which the critical current density is indepen- 
dent of H3 [5, 20]. At fixed H2 and variable H3, the 
slope of Sf(H 3 ) abruptly changes sign at the values of 
H2 , i?3 corresponding to the exact in-plane field orienta- 
tion. Similar results were obtained in the (Hi,H§) plane 
(not shown). This analysis gives the misalignment an- 
gles between the field components Hi and H2 and the 
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FIG. 2: Contour plot of amplified cavity transmission for 
variable H3 field while H2 = T (blue/gray: fiW, white: 
6.7 fiW). Scans maxima are shown by small squares, and are 
repeated for various H2 (the other contour plots are omitted 
for clarity). The maxima amongst a group of scans of given 
H2 is shown by large dots and indicate an in-plane field po- 
sition achieved by a specific (ff|,ifj) combination. Similar 
measurements are performed for the (Hi,Hs) field combina- 
tion (not shown). Insert: quality factors Qi,2 for (Hi,H^) 
(squares) and (if|,ff|) (dots) fields, respectively. 



film plane, 6 1 = 6.85° ± 0.02° and 9 2 = 3.42° ± 0.02°, 
respectively. The observed equal slopes of 5f(H^) for all 
accessible fields if 1.2 indicate that the effect of the per- 
pendicular field is decoupled from the shift caused by the 
in-plane field. Moreover, the quality factors for the fre- 
quency scans for different in-plane fields (H^ 2 ,H^) are 
essentially field independent, as shown in inset of Fig. 2. 
This result provides one more evidence that at (if* 2 , #3 ) 
no dissipation caused by possible penetration of vortices 
occurs up to IT, and the film is in the Meissner state. 

The main experimental evidence for NLME is shown in 
Fig. 3 in which the frequency shift is plotted as a function 
of the in-plane field H p . The (if* 2 , if*) maxima amongst 
scans of fixed Hi 2 exhibit noticeable shifts towards lower 
frequencies. For instance, for the data presented in 
Fig. 2, the decrease from the zero-field /o =19.213 GHz to 
19.174 GHz at Hi = 1 T gives a relative shift of -0.2%. 
Similarly, a relative shift of —0.58% is obtained for the 
scans in the (Hi, f/3) plane at H 1 = 1 T (last data point 
in Fig. 3(a)). The quadratic dependence of 5f on H p 
shown in Fig. 3(a) is consistent with Eq. (3). 

The NLME angular dependence is measured by rotat- 
ing a constant field H p within the sample plane, while 
monitoring continuously the resonance frequency f(H p ). 
Transmission scans obtained for an 180° rotation of a 
field H p =1 T are shown as a contour plot in Fig. 3(b) 
(the resonance peaks are identified by triangles). The fre- 
quency shifts at IT corresponding to Fig. 3(a) are shown 
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FIG. 3: (a) Relative frequency shift for in-plane field combi- 
nations (if*, if*) (squares) and (if J, if*) (dots, taken from 
Fig. 2) as a function of Hp (lines are through-zero fits), 
(b) Amplified transmission contour plot as a function of tp 
(blue/gray: fiW, white: 11 /iW) at constant in-plane field 
H p (ifi) =1.0 T. Maxima, shown by symbols, are in excellent 
agreement with the theory (black line). 



with their original symbols and indicated by the lateral 
arrows. The 180° periodicity predicted by Eq. (3) is ev- 
ident. The lines show the fit to the theory described 
below. 

To explain our data quantitatively we calculate A (if) 
in the dirty s-wave limit by solving the Usadel equa- 
tions [21] for the energy integrated Green functions, 
f u = e lB sin a, = cos a, 

D ( Q 2 \ 

— I V 2 a — sin 2a \ — lo sin a — A cos a, (4) 



ui>0 



A = 2ttjT > sin a, J = -2TreNDTQ 



^sin 2 a, (5) 

CJ>0 



where D is the electron diffusivity, 7 a pairing constant, 
lob the Debye frequency and N the density of states. For 
weak current pairbreaking, a = ao + 5a where sin ao = 
A/Vcj 2 + A 2 , and Sa(z) oc Q 2 depends only on z for the 
planar geometry shown in Fig. 1. Then Eq. (4) gives 

5a" - k 2 Ja = (Q 2 /2)sin2cv - (25A/D) cos a , (6) 



where fc 2 = 2Vw 2 + A 2 /L> and Q 2 (z) = 
2pqzsmif. In the case of no suppression of A at the 
surface, 5a(z) in Eq. (6) is a quadratic polynomial of z. 
Calculating coefficients in this polynomial, the ac cur- 
rent 5J(z) oc and the field-induced gap correction 
5 A oc Q 2 from Eqs. (5) and (6), we obtain A 2 (if) = 



-(/>odg/27ryUo J_ d / 2 &J(z)dz. For T = 0, this yields 
A 2 (ff P ) = A 



/ 2ir[i H p £, 



(C v d* + C&) 



(7) 



C v = (tt/96+ 1/18tt)(2-cos2<£>), (8) 
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where C 5 = tt 2 /32 + 7/24 + 1/6tt and £ 2 = D/A. 

The parameters in Eqs. (2) and (7) can be inferred 
from independent measurements, except the absolute po- 
sition of the origin ip = 0° relative to the coil axis 
i?2- Using the mean- free path extracted from the nor- 
mal state resistivity, £ = p^ 1 x 3.7 pOcm 2 = 1.6 nm 
[22] , we estimate the dirty limit penetration depth, Ao ~ 
A(^/^) 1/2 (T c0 /T c ) 1 / 2 = 226 nm, the coherence length 
£o = V&(T c o/T c y/ 2 = 9 nm, and k = A /£ - 25 for 
the clean- limit values A ~ £ ~40 nm [23]. Here the fac- 
tor T c q/T c accounts for the difference between our sam- 
ple critical temperature and T c0 = 9.2 K of pure Nb. We 
also used the nominal value of the geometrical inductance 
L g = 1295 pH for which the measured and designed zero 
field values of / agree very well, and set G = 1. 

The agreement of our data with Eq. (7) is quite good, 
as shown in Fig. 3. Here the origin ip — 0° was adjusted 
by 2.5° to account for the misalignment between the strip 
axis x and the coil axis Hi- The angular dependence 
shown in Fig. 3(b) is in excellent agreement with the 
theory. Also, the two quadratic cavity pulls shown in 
Fig. 3(a) lead to through-zero fitted slopes of df/ (fH 2 ) = 
-5.4 x 1CT 3 1/T 2 for ip = 90°, and -2 x 1CT 3 1/T 2 for 
(p = 0°, consistent with the corresponding theoretical 
values -5.7 x 1CT 3 and -2.2 x 10~ 3 . For d = 65 nm, the 
in-plane H c \ w 0.7 T is of order of the maximum field 
used in this study. We have also observed similar NLME 
behavior for a 40 nm film with H c \ = 1.28 T. 

Demonstration of the NLME in conventional high-K Nb 
films shows that our method is robust and readily usable 
to probe unconventional pairing symmetries in other ma- 
terials. In this case, field rotation in the plane of a thin 
film strip results in the orientational dependence of X(H) 
which, in addition to the geometrical factor 2 — cos2ip, 
would contain the intrinsic contribution from the orien- 
tational dependence of the order parameter. Separation 
of the two contributions will reveal the symmetry of the 
order parameter using the NLME angular spectroscopy 
proposed in this work. The observation of NLME on a 
superconductor with k = 25 reported here shows that 
our method is indeed insensitive to low bulk H c i values 
and can be applied to extreme type-II superconductors. 
The method can be adapted to study superconducting 
films made of complex materials, placed on a Nb cavity. 

Both types of frequency shifts shown in Figs. 2 and 3 
— quadratic NLME due to parallel field and the linear 
dissipative caused by perpendicular field need to be well 
understood in research areas requiring precise knowledge 
of the cavity pull, particularly in quantum computing 
studies. Quantum spins arranged in small crystals can be 
entrapped atop of on-chip structures [24] and show rel- 
evant multi-photon coherent phenomena [25], but often 
require strong magnetic fields to control their dynamics. 
The studies of on-chip superconducting cavities need to 
take into account the effects reported here. 

We have demonstrated the field-dependent non-linear 



Meissner effect under true equilibrium conditions in 
strong dc magnetic fields. The use of thin film resonator 
cavities opens up opportunities of probing intrinsic sym- 
metries of superconducting order parameter in supercon- 
ductors using NLME angular spectroscopy. 
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